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Abstract —The increasing penetration of renewable energy in 
recent years has led to more uncertainties in power systems. 
These uncertainties have to be accommodated by flexible re¬ 
sources (i.e. upward and downward generation reserves). In this 
paper, a novel concept. Uncertainty Marginal Price (UMP), is 
proposed to price both the uncertainty and reserve. At the same 
time, the energy is priced at Locational Marginal Price (LMP). A 
novel market clearing mechanism is proposed to credit the gener¬ 
ation and reserve and to charge the load and uncertainty within 
the Robust Unit Commitment (RUC) in the Day-ahead market. 
We derive the UMPs and LMPs in the robust optimization 
framework. UMP helps allocate the cost of generation reserves 
to uncertainty sources. We prove that the proposed market 
clearing mechanism leads to partial market equilibrium. We find 
that transmission reserves must be kept explicitly in addition to 
generation reserves for uncertainty accommodation. We prove 
that transmission reserves for ramping delivery may lead to 
Financial Transmission Right (FTR) underfunding in existing 
markets. The FTR underfunding can be covered by congestion 
fund collected from uncertainty payment in the proposed market 
clearing mechanism. Simulations on a six-bns system and the 
IEEE 118-bus system are performed to illustrate the new concepts 
and the market clearing mechanism. 

Index Terms —Uncertainty Marginal Price, Cost Causation, 
Robnst Unit Commitment, Einancial Transmission Right, Gen¬ 
eration Reserve, Transmission Reserve 


Nomenclature 


Indices 

i,l,t indices for generators, lines, and time intervals 
TO, n index for buses 
rrii index of bus where unit i is located 

k index of the worst point for uncertainty 

Functions and sets 

symbol for the optimal value of a variable 

E feasible set for UC and dispatch 

U uncertainty set 

Cf (•), C/(•) cost related to dispatch and UC for unit i 
£{■) Lagrangian function 

G{m) set of units located at bus to 

JC set of the indices for 


of indices k for upward and downward 
UMPs at bus TO time t 

Constants 

Nd, Nt number of buses and time intervals 
dm,t aggregated equivalent load 

Fi transmission line flow limit 

F; m shift factor for line I with respect to bus to 
pmm pmax minimum and maximum generation outputs 
r“,rf ramping-up/down limits between sequential in¬ 

tervals 

i?“, Rf ramping-up/down limits for uncertainty accom¬ 
modation 

Um,t bound for uncertainty 

e is the k*^ worst uncertainty vector in 1C, C 

e M 

FTRm_>„ FTR amount from bus m to n 
Variables 

/i t unit on/off status indicators 

unit start-up and shut-down indicators 
Pi^t generation dispatch 

P”^t net power injection 

em,t uncertainty at bus to time t 
Z optimal value of problem (SP) given 

{x,y,z,I, P) 

APi^t generation re-dispatch 

Aff°^ transmission capacity reserve in positive direction 

A/“® transmission capacity reserve in negative direc¬ 
tion 

AP^\ net power injection change 

Xt, X^ Lagrangian multipliers 

a,P,r] non-negative Lagrangian multipliers 

TTm t marginal prices, tt^ ^ for energy price; is the 

UMP for fcth uncertainty point; for upward 
UMP; for downward UMP 

QTt 1 Oi ““"upward and downward generation reserves 
T*™ t charge for uncertainty source 

Qff, credits to generation reserve for unit i and trans¬ 

mission reserve for line I at time t 
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I. Introduction 

N modern power systems, uncertainties grow significantly 
with the increasing penetration of Renewable Energy 
Source (RES), such as wind power generation. They pose 
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new challenges for the operation of electricity markets. In the 
Day-ahead market (DAM), the Unit Commitment (UC) and 
Economic Dispatch (ED) problems considering uncertainties 
have become a focus of research in recent years. The objective 
of the UC problem is to find the least cost UC solution 
for the second day while respecting both system-wide and 
unit-wise constraints. By fixing the UC variables, the ED 
problem is established. The Locational Marginal Price (LMP) 
and reserve price are then obtained as byproducts of the 
ED problem [1], [2]. When considering the uncertainties, the 
generation from uncontrollable RES are uncertain parameters 
in the optimization problem. 

Recently, Robust UC (RUC) is proposed to address the 
issues of uncertainty [3]-[7]. The largest merit is that the 
UC solution can be immunized against all the uncertainties 
in predefined set. The key idea of the two-stage RUC is to 
determine the optimal UC in the first stage which leads to 
the least cost for the worst scenario in the second stage. 
However, this approach is conservative and the Robust ED 
(RED) is absent. Authors in [8] combined the stochastic and 
robust approach using a weight factor in the objective function 
to address the conservativeness issue. [9], [10] employed the 
Affine Policy (AP) to formulate and solve the RED problem. 
A Multi-stage RUC is proposed to incorporate the latest 
information in each stage [11], where AP is also used to 
overcome the computational challenge. Recently, we reported 
a new approach which tries to bridge the gap of RUC and 
RED [7], [12]. 

In DAM, the main difficulty for pricing is that RED is 
absent in the traditional RUC [4], [5], [8]. On the other hand, 
a large number of works on pricing reserves exists within 
the UC framework considering contingencies [2], [13], [14] 
and stochastic security [15]. They are normally modeled as 
co-optimization problem. In [2], the reserve is cleared on 
zonal levels. Instead of countable contingency scenarios [15] 
or single additional scenario for reserve [16], the infinite 
continuous uncertainties are considered in the RUC, and the 
reserves are fully deliverable in infinite scenarios. In this 
paper, we propose a novel mechanism to price the energy, 
uncertainty, and flexibility simultaneously based on the RUC 
in [7]. An explicit price signal is derived for pricing the 
uncertainty. As the ED solution obtained is robust [7], both 
marginal impacts of the uncertainty and flexibility are reflected 
in these prices. In the proposed mechanism, reserve costs 
are allocated to uncertainty sources. Generation reserves, also 
called flexibilities in this paper, are the key factor for the 
robust optimization approaches. They are entitled to proper 
credits based on their contribution to uncertainty management. 
According to the market equilibrium analysis, market partic¬ 
ipants (price takers) can get the maximal profit by following 
the ISO/RTO’s dispatch instruction. 

The generation reserve and its deliverability are the main 
focus in [7]. The definition of LMP in [17] are employed 
to derive the energy price. The new concept. Uncertainty 
Marginal Price (UMP), is proposed to define the marginal cost 
of immunizing the next increment of uncertainty at a specific 
location. Load and generation are a pair, and they are priced 
at LMP. Uncertainty and flexibility (i.e., generation reserve) 


are another pair, and they are priced at UMP. Both LMPs and 
UMPs may vary with the locations due to transmission con¬ 
gestions. Limited by the transmission capacity and power flow 
equations, sometimes the uncertainties at certain buses cannot 
be mitigated by the system-wide cheapest generation reserve, 
and expensive generation reserve, which is deliverable, has 
to be kept in the system. Therefore, uncertainty sources are 
charged and generation reserves are credited based on UMPs 
at the corresponding buses. 

As the transmission reserve is kept within the RUC frame¬ 
work, the congestion component may exist in both the energy 
price and reserve price even if the physical limit of the line 
is not reached yet in the base case scenario. LMP congestion 
costs are distributed to Einancial Transmission Right (ETR) 
holders in the existing market according to the LMP differ¬ 
ence and the ETR amount. The revenue inadequacy occurs 
when the LMP congestion cost collected is smaller than the 
credit distributed to ETR holders, which is also called ETR 
underfunding. This has been a serious issue in recent years 
in the industry [18], [19]. We reveal that transmission reserve 
will be another reason for ETR underfunding when physical 
transmission limit is adopted in Simultaneous Eeasibility Test 
(SET) for ETR market [18], [20], [21]. This conclusion is 
applicable to any robust UC framework for DAM. 

The main contributions of this paper are listed as follows. 

1) The novel UMP for uncertainties and generation re¬ 
serves, as well as LMP for energy, are derived within a 
robust UC framework. The derivation is for uncertainties 
set with interval and budget constraints. The general 
concepts still apply when other uncertainty sets are 
modeled. 

2) It is revealed that transmission capacities have to be 
reserved for uncertainty accommodation and the trans¬ 
mission reserves may cause ETR underfunding because 
of the deficiency of energy congestion revenues based 
on existing market rules. 

3) A new market clearing mechanism is proposed to credit 
the generation and reserve and to charge the load and 
uncertainty. The payment collected from uncertainty 
sources can exactly cover the credits to generation 
reserves and transmission reserves, effectively resolving 
the ETR underfunding issue. 

The rest of this paper is organized as follows. Derivation of 
the LMP and UMP is presented in Section II, so is the market 
clearing mechanism for charge and credit based on LMP and 
UMP. Case studies are presented in Section III. Section IV 
concludes this paper. 

II. RUC AND Market Clearing 

One motivation of this work is to price the uncertainty, 
and allocate the cost of uncertainty accommodation to the 
uncertainty source. As the uncertainty source is charged the 
uncertainty payment, it has the incentive to reduce the uncer¬ 
tainty. With UMP, we can follow the cost causation principle, 
which is normally required in the market design, to charge the 
uncertainty sources. Cost causation principle is described as 
“require that all approved rates reflect to some degree the costs 
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actually caused by the customer who must pay them” in KN 
Energy, Inc. V. FERC, 968 F.2d 1295, 1300 (D.C. Cir. 1992). 

Another important motivation is to provide a theory that 
supports the application of the RUC in the DAM clearing. 
Although the RUC/RED are studied extensively, the only 
application of the RUC now is for the Reliability Assessment 
Commitment (RAC) in the DAM. There are several reasons 
why they are not applied in the DAM clearing. First, the 
computation burden of RUC is much larger than the standard 
UC. Second, as the objective is the cost of the worst-case sce¬ 
nario [3], [5], the solution is criticized on over conservatism. 
Third, no economic dispatch and prices are available within 
the RUC framework. Recently, with the new achievements 
in the algorithms, models, and high-performance computing 
application [6]-[8], [22]-[24], the first two obstacles are being 
addressed with great promises. This paper tries to clear the last 
obstacle with the new model [7]. Adopting RUC in the market 
clearing can give clear price signals for the uncertainties and 
reserves. On the other side, it is also easier for the solution to 
pass the robustness test, which is a RUC, in RAC. To our best 
knowledge, this is the first work on pricing energy, uncertain¬ 
ties, and reserves within the robust optimization framework in 
DAM. Hence, we focus on illustrating the concept with the 
following assumptions. 

• Network loss is ignored. Shift factor matrix is constant. 

• Uncertainty is from load and RES. Contingency is ig¬ 
nored. 

• The uncertainty budget set can be truly formulated by the 
ISO/RTO. 

A. RUC and RED 

ISOs/RTOs desire to get the optimal UC and ED solution 
in the base-case scenario. They can re-dispatch the flexible re¬ 
sources, such as adjustable load demands and generators with 
fast ramping capabilities, to follow the load when deviation 
occurs (or uncertainty is revealed). Consistent with the robust 
literature [4], [5], the uncertainty set is modeled as 

U := {ee ^ 

^ < Af,Vf} 

Af‘ is the budget parameter and assumed as an integer [3]. It is 
noted that all the flexible resources are modeled as generators. 
In this paper, the RUC is formulated according to the model 
in [7]. 

(RUC) min C^{x) + C^{p) 

(x,p)^T 

s.t. Ax Bp < b (1) 

J- := I (tc, p) : Ve G 3Ap such that 

Cx + Dp + GAp<d + Eey (2) 

The basic idea of the above model is to find a robust UC and 
ED for the base-case scenario. The UC x and dispatch p are 
immunized against any uncertainty e gU. When uncertainty e 
occurs, it is accommodated by the generation adjustment Ap. 
Please refer to Appendix A for the detailed formulation. 


Column and Constraint Generation (CCG) based method is 
used to solve the above model [6]. Problem (MP) and (SP) 
are established as follows. 

(MP) min C''^(a:) 3-C'^(p) 

{x,p) 

s.t. Ax + Bp < b 

Cx + Dp+ GAp'^ <d + Ee’^yk G JC (3a) 

and 


(SP) 

Z := max min l”'~s 

(4a) 


e^U (s, Ap)G7?.(e) 


7^(e) 

:= |(s, Ap) : s > 0 

(4b) 


GAp — s < d — Gx — Dp + Ee^ 

(4c) 


where K is the index set for uncertainty points e which are 
dynamically generated in (SP) with iterations. Please refer to 
Appendix B for the detailed formulation. It should be noted 
that is the extreme point of U. Variable Ap^ is associated 
with e^. The objective function in (SP) is to find the worst 
point in U given {x,p). The procedure is 
1: /C ^ 0, fc ^ 1, Z <—hoo, define feasibility tolerance <5 
2 : while Z > 5 do 

3: Solve (MP), obtain optimal {x,p). 

4: Solve (SP) with x — x,p — p, get solution (Z, e^) 

5: /C3—/CUfc, /c-l-l 

6 : end while 

Once the procedure is converged, we also get the optimal 
UC and ED solution by solving (MP). Similar to traditional 
LMP calculation, we fix the binary variables as x. Then a 
convex linear programming problem (RED) can be formed as 


(red; 

) min 
P,AP 

t i 

(5) 

s.t. 

(At) 


,t ~ ^ ^ 

(6a) 

(A.t) 

W,t < 

i,,tPr\yi,t 

(6b) 

(A.t) 

-p^,t 

< -WtPrwiA 

(6c) 

) 

p^,t- 

Pt,t-i < u“(i - yi,t) + pr''w,t,yi,t 

(6d) 

) 

-Pr,t 

+ Pt,t-i < rf{l - Zi^t) + PW^Zi^tWiA 

(6e) 

im,t) 


,^P'^,t<FiWlA 

(6f) 

{m,t) 

-E 

Ti,^P2t<FiWlA 

(6g) 

(Af) 

Y,APlt = Y.Wn,tWtWkGJC 

(7a) 

Wit) 

Pi,t + 

APit < WtPr^Whtwk G 1C 

(7b) 

Wit) 

-p^,t 

-APl^ < -WtPrWiAWkGlC 

(7c) 

Wit) 

^Plt 

<Rni-W,t)WhtWkelC 

(7d) 

Wit) 

-APi 

'^^t<Rfil-Zi,t+i),yt,tWkelC 

(7e) 

Wit) 


^miPZ + ^p'Wf)<Fi,yi,t,ykGic 

(7f) 

Wit) 

-E 

Ti^^iP^ + ap^A) < Fi,yi,t,yk G /C(7g) 


m 
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where ( 6 a)-( 6 g) are the constraints for the base ED, and (7a)- 
(7g) are constraints for different extreme points in U. (7a) 
denotes the load balance after re-dispatch. The generation 
adjustments respects capacity limits (7b)(7c) and ramping 
limits (7d)(7e). Network constraints are denoted by (7f)(7g). 
The and are defined as 

■Pm,*— X] Pi,t- 

ieO{m) 

and 

ieQim) 


respectively. 


B. Marginal Prices 


In this section, marginal prices for the energy, uncer¬ 
tainty, and generation reserve are derived based on the La- 
grangian function. Denote the Lagrangian function for (RED) 
as C{P, AP,X,a, P,ri), which is shown in Appendix C. Ac¬ 
cording to the definition of marginal price [17], the LMP for 
energy at bus m is 


dC{P,AP,X,a,P,T]) 


dd 


( 8 ) 


'7n,t 


=Xf.- 


'y ] ^l,m (yi,t y '.y ] (vi^ 


— id 

t Ul,t 


I keK 


It is observed that the impact of the uncertainty is also reflected 
in the LMR 

The new concept, UMP for DAM, is defined as the marginal 
cost of immunizing the next unit increment of uncertainty. Eor 
e^, an extreme point of U, the UMP is 


TT 


u,fe 

m,t 


d£{P,AP,X,a,d,r]) 


A,^-Er/.m {nit-Bit) 

i 

(9) 


Both the uncertainty and generation reserve are priced at 
In the derivation of the worst point is the only 

concern. Therefore, the general principles in this paper still 
work when U is replaced with other sets. 

It should be noted that (9) is intermediate price signals. In 
order to get the aggregated UMPs, the following new sets are 
defined based on the sign of 


^m,* := {k ■■ > 0 } 


^down 


:= {fc : tt;;;" < 0 } ( 10 ) 


The aggregated upward and downward UMPs are defined as 


‘ m,t 


u,fc 


u.down 


u.fc 


( 11 ) 


respectively. In the following context, we will show how the 
aggregated UMPs are used. 


C. Market Clearing Mechanism 

With LMP and UMP, the charges and credits for the 
market participants become clear and fair in the DAM. Energy 
clearing is straightforward. The basic principle related to 
uncertainty and flexibility is that those who cause uncertainties 
(uncertainty sources), such as RES, pay based on UMP and 
those who contribute to the management of uncertainties 
(uncertainty mitigators), such as generators or storage with 
ramping capabilities, get paid. 

1) Energy Payment and Credit: LSEs pay based on the 
amount of the load and LMP. The energy payment from the 
LSE at Bus m at f is tt^ tdm,t- It should be noted that RES 
is entitled to the credit due to the negative load modeled in 
RUC. Generator i, located at Bus m*, is entitled to the credit 
^m tPi,t for energy production. 

2) Charge to Uncertainty Source: The uncertainty source 
can be charged as 

'^m,t = '^rn,Pm,t ( 12 ) 

k&K 


The uncertainty source pays based on the marginal price and 
the worst point e^. The uncertainty source is charged only 
when is non-zero, and it may have to pay more when 
the uncertainty becomes larger. The uncertainty point ^ can 
be upward (i.e. ^ > 0) or downward (i.e. ^ < 0). We 

have the following lemma regarding the relation between the 
signs of and 4 . 

Lemma 1. If ^ > 0, then > 0. If ^ < 0, then 

TTm.* < 0 . 


Please check Appendix D-A for the proof. When the budget 
set is adopted, the extreme point ^ € {—Um,t, 0, Um,t} [4], 
[7], so the uncertainty charge in (12) can also be written as 
(13) according to Lemma 1 and (11). 


iTf UjUp I u.aown/ > 

^7n,t — '^m “r t \ 


(13) 


Thus, upward and downward uncertainties are charged sepa¬ 
rately. It should be noted that we still need to use (12) when 
other uncertainty sets are used. 

3) Credit to Generation Reserve: Only resources that can 
provide deliverable generation reserve are entitled to credits. 
If * € G{m), then the credits can be formulated as 

= E <tA/^V (14) 

fce/c 


In other words, generation reserve is paid the UMP at the bus 
where it is located. If = 0, then the associated credit 
is zero no matter what the value of APf^^ is. Similar to the 
uncertainties, the generation reserves can be in either upward 
or downward direction. Denote the upward generation reserve 
as and the downward generation reserve as ( 3 ^°™ 

QZ := min - P.,*, i?“(l - y,,*)} , (15) 

QZr ■■= - Pm, -Pfil - 2 i,*+i)} .(16) 

We also have the following lemma regarding the relation 
between QZ, 0 *°™ ^Pit- 
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Lemma 2. If i € Gim), then the optimal solution to 

problem (RED) is 


and 





up 

down 
i,t ’ 


if TT. 

if TT. 


.U,fc 

m,t 

,u,k 

7n,t 


> 0 
< 0 




(17) 


Please check Appendix D-B for the proof. The credit to 
generation reserve i located at bus m (14) can be rewritten as 
(18) according to Lemma 2 and (11). 


r)G _ u.up^up , u.down^down 


(18) 


(18) shows that the upward and downward generation reserves 
are credited separately. Flexible resources may receive credits 
for both the upward and downward generation reserves simul¬ 
taneously. (18) always holds even if other uncertainty sets are 
modeled in RUC. 


D. Transmission Reserve and Revenue Adequacy 

Some transmission capacities are reserved according to the 
solution to (RED). These transmission reserves are used to 
ensure the ramping deliverability when the uncertainty is 
revealed, as shown in (7f) and (7g). It is noted that they 
are determined automatically in (RED), and kept explicitly 
without explicit transmission reserve requirement constraints. 
Just like the “scheduled” generation reserve, the “scheduled” 
transmission reserves in positive direction and negative direc¬ 
tion are 

■■=Fl-Y.Ti,raP2u ( 19 ) 

m 

( 20 ) 

m 

respectively. They are always non-negative. 

An important issue related to the transmission reserve is 
the credit entitled to the Financial Transmission Right (FTR) 
holders. FTR is a financial instrument used to hedge con¬ 
gestion cost in the electricity market, where participants are 
charged or credited due to the transmission congestion [21], 
[25]. Within the robust framework, the effective transmission 
capacity for base-case scenario is different from the physical 
limit, which is used in the Simultaneous Feasibility Test (SFT) 
for FTR market [18], [20], [21]. In the existing market, the 
FTR credit is funded by the energy congestion cost, which 
is the net payment of energy. However, the energy congestion 
cost may not be sufficient to fund the FTR credit [18], [19]. We 
argue that the transmission reserve becomes a new reason for 
FTR underfunding in any framework to guarantee the ramping 
deliverability. 

Theorem 1. If transmission reserve Aff°^ and are 

kept for line / at time t in DAM, then the maximum FTR 
underfunding associated with line I at time t is 

kGJC 

due to the dehciency of energy congestion cost. 



Fig. 1. Money flow of the proposed market clearing mechanism, where 
uncertainty sources make the uncertainty payment, and LSEs make the energy 
payment. 


Please check Appendix D-C for the proof. From the FTR 
holder’s point of view, (21) is the credit due to the transmission 
reserve. Therefore, we also call (21) transmission reserve 
credit, and denote it as 

( 22 ) 

keic 

At most one of fjf^ and is non-zero for transmission 1. 
The credit to positive transmission reserve is zero for line I at 
time t, when either J2k^ic ^it = 0 or Afi°^ = 0. 

Theorem 2. If (RED) is feasible, then uncertainty payment 
can exactly cover generation reserve credit and transmission 
reserve credit, and the revenue adequacy is always guaranteed 
in the proposed market clearing mechanism. 

Please check Appendix D-D for the proof. Theorem 1 
reveals that FTR underfunding issue can occur within the 
existing market structures as long as the transmission reserve 
is non-zero, even if the LMPs are calculated based on other 
approaches. Theorem 2 shows that the new market clear¬ 
ing mechanism overcomes the FTR underfunding issue. The 
money flow of the proposed market clearing mechanism is 
depicted in Fig.l. Energy payment collected based on LMP 
is distributed to FTR holders as LMP congestion cost and 
generators as energy credit. On the other hand, the payment 
collected based on UMP is distributed to FTR holders as 
transmission reserve credit and flexible resources as generation 
reserve credit. The LMP congestion cost and transmission 
reserve credit can exactly cover the FTR credit, which is 
calculated based on the LMP difference and FTR amount. 


E. Market Equilibrium 

In this section, we characterize the competitive market equi¬ 
librium model. In the electricity industry, the partial market 
equilibrium model is often employed [1], [13], [26], where 
market participants are price takers [27]. 

The energy is cleared according to (6a). Uncertainty and 
generation reserve are cleared according to (7a). Without loss 
of generality, consider unit i located at bus m. Its profit 
maximization problem can be formulated as 


(PMPi)max 




u,down /^( 

fry, + 


Up^Uf 

''m,t 

-cr{p.,t) 

s.t. (6b) - (6e), (15) - (16) 


down ' 
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where the decision variable is * given the price signal 
(tt^ As proved in Appendix D-E, unit i 

is not inclined to change its power output level as it can 
obtain the maximum profit by following the ISO’s dispatch 
instruction Price signal ^ provides the incentives for 
unit i to dispatch power output to Pi t, and price gives 
incentives for unit i to maintain the generation reserve for 
uncertainty. Hence, the dispatch instruction P^ * and price 
signal (tt^ constitute a competitive partial 

equilibrium [27]. 

F. Discussions 

As the Pi^t and (or {Q'^^t,Qf°™)) are coupled by 

(7b) and (7c), the opportunity cost — f3ft) is enough to 
provide the incentives for i to keep the generation level at Pi^f 
Including — in the generation reserve price has several 
benefits. Firstly, generation reserves provided by different units 
are priced fairly. Generation reserve prices are the same for the 
units at the same bus, and they may vary with locations if line 
congestions exist. Secondly, higher generation reserve price 
attracts long-term investment for flexible resources. Thirdly, it 
is consistent with the existing reserve pricing practice [ 2 ], [28]. 
In fact, generation reserve price is consistent with the UMP. 
Therefore, the uncertainties and flexibilities are also treated 
fairly at the same bus. 

The upward and downward UMPs are obtained according to 
(11), respectively. The uncertainty sources are charged accord¬ 
ing to (13). The generation reserves are credited according to 
(18). The price signal P^t defined in (9) and re-dispatch AP/^^ 
are intermediate variables for market clearing. The proposed 
UMP may be non-zero even if the uncertainty at a bus is zero. 
This is similar to the LMP, which may also be non-zero for 
the bus without load. 

The market clearing mechanism proposed in this paper 
follows the cost causation principle for the cost allocation. 
In reality, it may be controversial to allocate the reserve cost 
to uncertainty sources. However, we argue that it would be fair 
and must be done when the RES penetration level is high. An 
extreme case is when the loads are all supplied by RES. There 
has been study showing it is possible that the increasing RES 
penetration can cause higher system operation cost. This issue 
cannot be handled by the existing market clearing mechanism, 
in which loads pay for the additional system reserve that is 
required to accommodate the uncertainty from RES. In other 
words, loads are actually providing subsidies to RES. When 
the RES penetration level is low, the subsidies can help the 
growth of the RES. However, when the RES penetration level 
is high, these growing subsidies will cause serious fairness 
issue. On the other hand, with UMP as the stimulating price 
signals, RES will have the incentives to improve its forecast 
techniques and reduce its uncertainty. In the ideal case when 
its uncertainty approaches zero, RES will no longer pay. 

Following the existing practice, the UC variables are fixed 
during the marginal price derivation. Hence, the uplift issue, 
which exists in the real market, still remains in the proposed 
market clearing mechanism. Although the UC variables are 
fixed, the LMP and reserve price in the real market can provide 


effective signals for the long-term investment of generation 
and transmission as well as consumption strategy of electricity. 
Similarly, the uncertainty impact is not only reflected in UC, 
but also in the ED within the RUC model in this paper. 
Hence, the proposed LMP and UMP can also provide signals 
for the long-term investment of flexibilities (i.e. generation, 
transmission, and demand). 

The pricing for uncertainties proposed in this paper is not 
in conflict with the pricing for traditional reserves, which 
are mainly prepared for the contingencies. The traditional 
reserve prices can be derived in the framework by adding extra 
traditional reserve constraints, and the corresponding reserve 
costs can still be allocated to LSEs. 

It is observed that the credit in (14) is the sum of credits for 
all extreme points. That is because the related constraints may 
be binding for multiple extreme points, and the dual variables 
(shadow prices) for these constraints work together in the dual 
problem. The traditional price for energy and reserve also has 
similar form when multiple contingencies are modeled. 

Although only one scenario will happen in reality, we still 
need to consider the worst scenario defined in uncertainty set 
and keep enough reserves in DAM. That is because DAM is 
a financial market, and the LMP and UMP are the financially 
binding prices. This is similar to the existing market model 
considering contingencies. Even if the contingency seldom 
occurs, they are still modeled for market clearing, and the 
contingencies are reflected in LMP and reserve price. 

The issue of price multiplicity still exists in the proposed 
model [29] because problem (RED) is a linear programming 
(LP) problem. However, the price is unique with the nonde¬ 
generacy assumption. For simplicity, we have considered a 
single-sided auction in the proposed model. By introducing 
the demand bids, we can formulate a double-sided auction 
and the general principles in this paper will still apply. 

HI. Case Study 

A six-bus system and the IEEE 118-Bus system are sim¬ 
ulated to illustrate the proposed market clearing mechanism. 
In the six-bus system, the basic ideas of UMP are presented 
within the robust optimization framework. FTR underfunding 
issue is illustrated and a comparison between the UMP and 
traditional reserve price are presented. In the IEEE 118-Bus 
system, the UMP related products are presented for different 
uncertainty levels. The behaviors and impacts of flexible 
sources are analyzed by an energy storage example. 

A. Six-bus System 

A six-bus system is studied in this section. The one-line 
diagram is shown in Fig. 2. The unit data and line data 
are shown in Table I and Table II, respectively. Table III 
presents the load and uncertainty information. Column “Base 
Load” shows the hourly forecasted load. Assume that the load 
distributions are 20%, 40%, and 40% for Bus 3, Bus 4, and 
Bus 5, respectively. mi t and U 3 ^ in Table III are the bounds 
of the uncertainties at Bus 1 and Bus 3, respectively. The 
uncertainty bounds at other buses are 0 . 



7 



Fig. 2. One-line diagram for 6-bus system. 


TABLE IV 

Marginal Costs at Different Generation Levels ($/MWh) 



Gen, 

. 1 


Gen 

. 2 


Gen 

. 3 

TDW 

Dw 

mar. cost 

JDW 

i 2 

Dw 

^2 

mar. cost 

TDW 

^3 

pw 

mar. cost 

100 

124 

14.396 

10 

28 

32.638 

10 

12 

17.71 

124 

148 

14.588 

28 

46 

32.674 

12 

14 

17.73 

148 

172 

14.78 

46 

64 

32.71 

14 

16 

17.75 

172 

196 

14.972 

64 

82 

32.746 

16 

18 

17.77 

196 

220 

15.164 

82 

100 

32.782 

18 

20 

17.79 


TABLE I 

Unit Data for the 6-bus System 


# 

j^rnin 

^max 

Po 

a 

b 

c 


R'^ 

Cu 

Cd 

'j^on 

’J-’Off 

To 

1 

100 

220 

120 

0.004 

13.5 

176.9 

24 

24 

180 

50 

4 

4 

4 

2 

10 

100 

50 

0.001 

32.6 

129.9 

12 

12 

360 

40 

3 

2 

3 

6 

10 

20 

0 

0.005 

17.6 

137.4 

5 

5 

60 

0 

1 

1 

-2 


pmm pmax pp. jnin/max/initial generation level (MW); 
fuel cost ($): aP^ + bP + c ; 

ramping up/down rate (MW/h); 

Cu,Cd’ startup/shutdown cost ($); 

pon poff p^. on/min off/initial time (h) 


It is assumed that the relative forecasting errors increase 
with hours. Uncertainty ei^t and also respect 

A • ^ ^ A • j Ift (23a) 

y"l#^<AA,Vf, (23b) 

where (23a) denotes the uncertainty interval at a single bus, 
and (23b) represents the system-wide uncertainty [5]. The A 
and A"^ are the budget parameters for the single bus and 
system, respectively. 

1) LMP and UMP: Consider the case where A = 1, A'^ = 
2. The CCG based approach converges after 2 iterations. 


TABLE II 

Line Data eor the 6-bus System 


from 1 1 2 5 3 2 4 

to 2446635 

x(p.u.) 0.17 0.258 0.197 0.14 0.018 0.037 0.037 

capacity(MW) 200 100 100 100 100 200 200 


TABLE Ill 

Load and Uncertainty Data for the 6-bus System (MW) 


Time (h) 

Base Load 


U3,t 

Time (h) 

Base Load 

Ui,t 

y'3,t 

1 

175.19 

1.09 

0.29 

13 

242.18 

19.68 

5.25 

2 

165.15 

2.06 

0.55 

14 

243.6 

21.32 

5.68 

3 

158.67 

2.98 

0.79 

15 

248.86 

23.33 

6.22 

4 

154.73 

3.87 

1.03 

16 

255.79 

25.58 

6.82 

5 

155.06 

4.85 

1.29 

17 

256 

27.2 

7.25 

6 

160.48 

6.02 

1.6 

18 

246.74 

27.76 

7.4 

7 

173.39 

7.59 

2.02 

19 

245.97 

29.21 

7.79 

8 

177.6 

8.88 

2.37 

20 

237.35 

29.67 

7.91 

9 

186.81 

10.51 

2.8 

21 

237.31 

31.15 

8.31 

10 

206.96 

12.94 

3.45 

22 

232.67 

31.99 

8.53 

11 

228.61 

15.72 

4.19 

23 

195.93 

28.16 

7.51 

12 

236.1 

17.71 

4.72 

24 

195.6 

29.34 

7.82 


TABLE V 

Generation and Reserve (A = 1, A^ = 2, MW) 


T 

Pi 

P2 

Pa 

Q7 

Qfw" Q“P 

^down ^down 

^2 ^3 ^3 

21 

195.19 

25.58 

16.54 

24 

-24 12 

-12 3.46 -5 


Hence, JC = {1,2}. Given the UC solutions, the problem 
(RED) can be solved by commercial LP solver. The marginal 
prices are then obtained as byproducts. 

The generation outputs are presented in Table V at Hours 21. 
It can be observed that G1 supplies most of the loads at Hour 
21, which is 195.19 MW. According to the bid information in 
Table IV, G2 is much more expensive than G1 and G3. Hence, 
the output of G2 is relatively small and at the low level of 
its capacity. The upward and downward generation reserves 
provided by the three units are also listed in Table V. These 
data can be obtained directly from Eqs. (15) and (16) given 
the generation output (. Although the remaining generation 
capacity of G1 is 220 — 195.19 = 24.62 MW, the upward 
reserve is limited by its upward ramping rate 24 MW. In the 
meantime, the upward reserve provided by G3 is limited by its 
generation capacity although it has more remaining ramping 
capacity (i.e. min{20 — 16.54, 5} = 3.46 MW). 

Table VI shows the extreme points obtained in the CCG- 
based approach. The intermediate price signals for these points 
are also presented. It can be observed that the worst point 
is always obtained at the extreme point of the uncertainty set. 
Eor example, at Hour 21, the e} is 31.15 MW. It is exactly the 
upper bound of the uncertainty at Hour 21 at Bus 1. The data 
in Table VI also verifies Lemma 1. The intermediate UMPs 
have the same sign as the uncertainties j at the same 

bus. 

The LMPs, aggregated upward UMPs, and aggregated 
downward UMPs at Hour 21 are shown in Table VII. It is 
noted that UMPs still exist at buses without uncertainties (i.e.. 
Buses 2,4,5,6). This is similar to LMPs, which also exist 
at buses where net power injections are 0. The LMPs vary 
with locations, which indicates that the line congestion exists. 


TABLE VI 

Extreme Points of Uncertainty Set 



k 

= 1 



k = 

2 


J. 

^1 

U,1 

U,1 

"2 

^2 

u,2 

u,2 

* U.t 


^l,t 

^3,t 

U.t 



^3,t 

21 31.15 

8.31 

14.87 

14.87 

-31.15 

8.31 

-17.67 

1.77 



TABLE VII 

LMP AND UMP AT Hour 21 (A = 1, = 2) 


Price 

Busl 

Bus2 

Bus3 

Bus4 

Bus5 

Bus6 

TT^ 

14.97 

32.64 

34.4 

43.71 

41.94 

35.26 

7r“’“P 

14.87 

14.87 

16.63 

25.94 

24.17 

17.49 

.^Ujdown 

-17.67 

0 

0 

0 

0 

0 


The load at Bus 4 has to pay the highest LMP $43.71/MWh. 
The UMPs are also different at various locations. The highest 
upward UMP at Hour 21 is also located at Bus 4. With these 
prices, the market participants can be paid and credited. 

The LMP paid to G3 is $35.26/MWh on Bus 6, which 
is $17.49/MWh larger than its marginal cost. In the same 
time, The upward UMP is $17.49/MWh on Bus 6, which is 
exactly the difference between the LMP and G3’s marginal 
cost. Hence, G3 is the UMP setter on Bus 6. The UMPs 
provide important price signals on the planning of renewable 
energy sources and storages. For example, the UMP at Bus 2 is 
relatively small, so it is an ideal location for renewable energy 
sources in terms of payment for uncertainties. In contrast, 
the UMP at Bus 4 is large, which may attract the long-term 
investment for storages or generation plants with large ramping 
rates. 

2) Comparison with Existing LMPs and Reserve Prices: 
The motivation of this part is to compare the proposed clearing 
scheme with the existing one. However, as the reserve is not 
robust in the traditional scheme, we cannot compare them 
fairly. With the observation that the transmission constraints 
are the most challenging one in the robust UC framework, 
we drop these constraints in this subsection and add reserve 
constraints as follows. 

h,tPr < Qtr + QZ + p^,t < htPr, vi, t (24a) 
-Rfh,tAT < QZt < R^h.t^T, Vi, t (24b) 

Y,Qtr<Ru (24c) 

i i 

where QZt and are the largest upward and downward 

reserves, respectively. Rt and Rt are system-wide reserve 
requirements. Refer to [2], [30] for more details on the 
reserve formulations. In the experiment, AT is set to 1 and 
A = 0.8, = 2. The reserve requirements Rt and Rt are set 

to the lower and upper bounds of the system-wide uncertainty 
in (23b), respectively. 

The results are as expected. The optimal solutions of the 
RUC and the standard UC with explicit reserve constraints are 
the same. LMPs calculated in the RUC and UC also have the 
same values. The UMPs calculated in the proposed mechanism 
are also exactly the same as the reserve prices in standard 
UC. Two things are verified with these results. First, without 
transmission constraints, the solution to standard UC can easily 
be robust by adding reserve constraints. Second, the proposed 
LMPs and UMPs are consistent with LMPs and reserve prices 
in the existing market when the transmission constraints are 
dropped. 

When considering transmission constraints, the generation 
reserve cannot be guaranteed at bus levels in the traditional 




Bus 


Bu.s 


(a) Hour 20 (b) Hour 21 

Fig. 3. Upward UMP (blue bar) and reserve price (red bar) with network 
constraint. 


UC model. For simplicity, we assume that the 6 buses are 
in a zone. Consider the case where A = 0.8, A^ = 2. The 
upward UMP and reserve price at Hour 20 are depicted in 
Fig. 3a. It is observed UMPs at Bus 1 and Bus 2 are lower 
than the traditional reserve prices. In the same time, the UMPs 
at Bus 4 and Bus 5 is higher than the traditional reserve 
prices. The differences are caused by the congestion of Line 
1-4 for reserve delivery. It is worth mentioning that the LMP 
differences in two models are within 1% at Hour 20. The 
prices illustrated in Fig. 3b reveals another trends that the 
UMP may be higher than the traditional reserve prices. At 
Hour 21, the zonal reserve price is 0 while the UMPs are non¬ 
zeros at bus 4, 5, and 6. Because the constraint related with 
reserves in the RUC is stronger than the one in traditional UC 
model. Consequently, more expensive resources are used in 
RUC, which also generally leads to higher UMPs. 

3) FTP Underfunding: When A'^ = 2, A = 1, the gen¬ 
eration schedules at Hour 21 are 195.193MW, 25.577MW, 
and 16.54MW. The power flow of Line 2 is 97.63MW, 
which is 2.47MW smaller than its physical limit of lOOMW. 
The transmission reserve 2.47 MW is kept to guarantee 
the delivery of the generation reserve. The binding con¬ 
straint for Line 2 causes LMP differences. Hence, the 
FTR holder gets credits. Consider a set of FTR amounts 
[202.3429,23.2771,-55.772,-94.924,-94.924,20]. It can 
be verified that the FTR amounts satisfy the SFT in the FTR 
market. Then the total credit for the FTR holders is $5,554.77. 
However, the congestion cost in the DAM is $5,422.87. It 
means that the LMP congestion cost collected is not enough 
to cover the FTR credit. The FTR underfunding value is 
$5554.77 - $5422.87 = $131.90. 

The revenue residues after UMP settlement is $131.9. It ex¬ 
actly covers the FTR underfunding in this scenario. Therefore, 
the revenue is adequate at Hour 21. 

B. IEEE 118-Bus System 

The simulations are performed for the IEEE 118-bus system 
with 54 thermal units and 186 branches in this section. The 
peak load is 6600MW. The detailed data including generator 
parameters, line reactance and ratings, and load profiles can be 
found at http://motor.ece.iit.edu/Data/RUC118UMP.xls. Two 
cases are studied in this section. 

1) The uncertainty levels and load levels are changed to 
analyze the simulation results in the system level. The 
impact of transmission line capacity on prices is also 
studied. 
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TABLE VIII 

Operation Cost and UMP Payment ($,A'^ = 10) 


A 

Op. Cost 

Un. Payment 

Gen. Res. Credit 

Rev. Res. 

0.2 

1,866,023 

11,043 

10,560 

483 

0.25 

1,871,364 

20,044 

19,209 

835 

0.3 

1,877,471 

30,879 

28,658 

2,221 


2) An energy storage is installed at a specified bus with 
high UMP to show the potential application of UMPs. 

1) Case 1: We assume that the uncertainty sources are lo¬ 
cated at buses (11,15,49, 54, 56, 59,60,62,80,90). The bud¬ 
get parameter is set to 10 in this section. The bus- 
level uncertainty budget parameter A changes from 0.2 to 
0.3, and the bound of the uncertainty is the base load. The 
simulation results are shown in Table VIII. It can be observed 
that the total operation cost increases with increasing A. It 
indicates that a larger uncertainty level may increase the 
operation cost. The columns “Un. Payment” and “Gen. Res. 
Credit” denote the total payment from uncertainty sources and 
credit to generation reserves, respectively. The lowest payment 
is $11,043 and the highest one is $30,879. On the other 
hand, the credit entitled to the generation reserves is also a 
monotonically increasing function of A. When A = 0.3, the 
generation reserves have the highest credit. The last column 
“Rev. Res.” shows that the revenue residues related to UMPs. 
It can be observed that the residue is always positive. 

Fig. 4a in the next page depicts the heat map for the 
upward UMPs from Bus 80 to Bus 100 in 24 hours. The x- 
axis represents time intervals and the y-axis represents bus 
numbers. The color bar on the right shows different colors for 
various UMP values. For example, the $0/MWh is denoted by 
the blue color at the bottom, and the $18/MWh is represented 
by the dark red color on the top of the color bar. It can be 
observed that the uncertainty sources have system-wide unique 
UMPs at some intervals, such as Hours 8, 13, 15, and so on. It 
indicates that there is no transmission reserve in these hours. 
On the other hand, the UMPs at Hour 11 vary dramatically 
with different locations. The highest upward UMP is around 
$18/MWh, and the lowest one is around $2/MWh. According 
to the data shown in Fig. 4a, the high UMP at Bus 94 
may attract investment of flexible resources, such as energy 
storages, in terms of generation reserve credit, and Bus 100 is 
an attractive location for the investment of renewable energy 
sources in terms of uncertainty payments. 

Fig. 5 shows the uncertainty payment and generation reserve 
credit with respect to load levels. The base load level is set 
at 100%. Higher loads in general lead to more uncertainty 
payments and generation reserve credits. It is also consistent 
with the heat map of UMPs in Fig. 4a, where UMPs at peak 
load hours are high. It suggests that the generation reserves 
also become scarce resources when load levels are high. 

The transmission line capacity plays an important role in the 
price calculation. Fig. 6 shows the UMPs and upward UMPs 
at Hour 11 with respect to increasing capacity of Line 94- 
100. The prices at Buses 88, 94, and 100 are depicted. When 
the line capacity increases from 165MW to 175MW, LMP at 



Load Level (%) 

Fig. 5. Uncertainty payment (UP), generation reserve Credit (GRC), and 
operation cost (OC) with different load levels (A = 0.25) 


L]y[p Upward UMP 



Fig. 6. LMP (left) and upward UMP (right) at Hour 11 with respect to 
increasing capacity of Line 94-100 

Bus 94 decreases from $47.92/MWh to $35.84/MWh and that 
at Bus 88 also drops to $30.58/MWh from $38.78/MWh. The 
upward UMPs at Bus 94 and Bus 88 also drop by $8.20/MWh 
and $12.08/MWh, respectively. In contrast, the LMP and 
upward UMP at Bus 100, which is connected to Line 94-100, 
remain at $19.42/MWh and $L64/MWh, respectively. It shows 
that the change of line capacity may only have impacts on the 
prices at some buses. When the line capacity further increases 
to 185MW from 175MW, the changes of LMPs and UMPs at 
Bus 94 and Bus 88 are within $0.1/MWh, and there is still no 
change at Bus 100. It means that the additional lOMW cannot 
help deliver cheaper energy and reserves to Bus 94 and Bus 
88 . These results are also consistent with the analysis of the 
traditional LMPs [31]. 

2) Case 2: As discussed in Case 1, the upward UMP on 
Bus 94 is high at Hour 11. Assume that an energy storage 
(8MW/30MWh) is installed at Bus 94. A simple model for 
the energy storage is formulated as follows. 

0 < Pt < 

0 < -Pt < 

0 < pf < 

+ /f <l,Vf 

PNt Pq : 

where Et denotes the energy level, P/^ and Pf represent 
the discharging and charging rates, and ip and Ip are the 
indicators of discharging and charging. As the UMP is the 
major concern in this section, we use simplihed parameters for 
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(a) Without Storage at Bus 94 (b) With Storage at Bus 94 

Fig. 4. Heat Map for Upward UMPs (A = 0.3). Different colors represent various UMPs. Figure (a) depicts the UMPs from Bus 80 to Bus 100 in 24 hours 
without the energy storage at Bus 94. Figure (b) depicts the new UMPs after the energy storage is sited at Bus 94. 


storage. The discharging efficiency p‘^ and charging efficiency 
are set to 100%. The capacity and initial energy 

level i?o are set to 30 MWh and 15 MWh, respectively. The 
maximal charging rate and discharging rate are set to 
8 MW/h. 

By siting the energy storage, we can lower the new operation 
cost to $1,875,211 from $1,877,471. The payment collected 
from the uncertainty sources becomes $27,473, and the credit 
to generation reserves decreases to $24,289. Compared to the 
data in Table VIII, the energy storage also helps to reduce the 
payment related to UMPs. The storage is entitled to $1326 
generation reserve credit. Fig. 4b depicts the new upward 
UMPs after the installation of the energy storage. Compared 
to that in Fig. 4a, the upward UMP for Hour 11 at Bus 94 
decreases a lot. The UMPs for Hour 10 and 12 are also lower. 
It suggests that sitting the energy storage at Bus 94 effectively 
lower the generation reserve price. 

The simulation results demonstrate that flexible resources 
can lower the UMPs, and UMPs provide the investment signal 
at locations where generation reserves are scarce resources. 

IV. Conclusions 

A novel market model in this paper clears uncertainty, 
energy, and generation reserve simultaneously within the RUC 
framework in DAM. The uncertainty sources are charged 
and the generator reserve providers are credited based on 
the proposed UMP. The UMP formulation is derived within 
a robust optimization framework. We also characterize the 
market equilibrium for the new market clearing mechanism. 
As the market clearing mechanism is established within the 
robust optimization framework, the robustness of the dispatch 
is guaranteed. The optimal reserves for uncertainty accommo¬ 
dation are obtained in the model. The UMP proposed in this 
paper can effectively address the issue on how to charge and 
credit the uncertainties and generation reserve fairly in the 
market with RES. 

Our study also shows that traditional pricing mechanism 
within RUC framework may lead to FTR underfunding. The 
proposed market clearing mechanism can address this issue. 
Our study shows load serving entities can have lower energy 


prices within the new market scheme, as the reserve fees are 
paid by uncertainty sources. 

Many potential applications on UMPs are open. As UMPs 
are unified prices of uncertainties and reserves, it is interesting 
to investigate the optimal strategy for the one who is the 
uncertainty source as well as reserve provider in the market 
(e.g. the wind generation company with energy storages). The 
UMPs derived in this paper also provide an important price 
signal for the long-term investment of flexible resources. When 
the upward UMP or downward UMP at a bus is high, the 
investor can get more return in terms of generation reserves. 

Another potential future research on UMP is to study how to 
determine the budget uncertainty set in the market. Modeling 
the traditional spinning reserve for the contingency [13], [30] 
is also our future work. In this paper, the demand bids are 
not considered. We have forecasted load, forecasted RES, and 
uncertainty of load and RES for market clearing with a single¬ 
sided model. In an extended double-sided model, we can have 
demand bids, forecasted RES, and uncertainty of load and 
RES for market clearing. The forecasted load, forecasted RES, 
uncertainty of load and RES can be used in RAC. 


Appendix A 

Detailed Eormulation for Problem (RUC) 


(RUC) min E (25a) 

t i 

(25b) 


s-t. 


E Pi,t - dm,t\ < Fi,\/l,t(25c) 

h,tPr < P^,t < h,tPr\ Vf, t (26a) 

Pi,t - Pi,{t-1) < u“(l - Vt,t) + Pr''yi,t,yi,t (26b) 
-Pt,t + Pi,(t- 1 ) < - Zi,t) + (26c) 

minimum on/off time limit 
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and 

T :=|(a:, j/, 2 , /, P) : Ve € U, 3AP such that 

APi,t = Em.t, vt, (27a) 

i m 

< P,.t + AP,,t < /,,tP”“, Vi, i (27b) 

-Pf(l - z,.t+i) < AP,,i < P”(l - y.,i),Vi,i (27c) 
E APi,t - Cm,(27d) 

ieG(m) 

-Fi < Y, ^i,^{PZt + AP”V) < Pi,, VZ, i}. (27e) 

m 

The basic idea of the above model is to find a robust UC and 
dispatch for the base-case scenario. In the base-case scenario, 
(25b) denotes the load balance constraint; (25c) represents the 
transmission line constraint; (26a) denotes the unit capacity 
limit constraint; (26b)-(26c) denote the unit ramping up/down 
limits. Ji t, j/i t, and Zi^t are the indicators of the unit being 
on, started-up, and shutdown, respectively. Units also respect 
the minimum on/off time constraints which are related to 
these binary variables [1]. The UC and dispatch solution are 
immunized against any uncertainty e gU. When uncertainty e 
occurs, it is accommodated by the generation adjustment AP^ t 
(27a). Generation dispatch is also enforced by the capacity 
limits (27b). (27c) models the ramping rate limits of generation 
adjustment AP^ (. In fact, the right and left hand sides of (27c) 
can correspond to a response time AT, which is similar to the 
10-min or 30-min reserves in the literatures [30]. (27e) stands 
for the network constraint after uncertainty accommodation. 

Appendix B 

Detailed Formulation for Problem (MP) and (SP) 


EE(^*''(^M)+c'/fc) 


(MP) min 

{x,y,z,I,P,AP) 

S.T. (25b), (25c), (26a)-(26c), minimum on/off time limit 


Y^Pt,t = Y^'rn,t,yt,ykGlC 


(28a) 


h,tPr < p^,t + ^pit < ktPryh c vfc e /c ( 28 b) 

AP^, <P“(l-j/,,t),Vz,f,Vfce/C (28c) 

< Rf{l - G JC (28d) 

-Fi < + AP”jf) < P,,Vfc G K,yi,t (28e) 

m 


A pinj 


,fc 


y] Ap,-et.z,Vm,f,VfcG/C, 

iep(m) 


(28f) 


^PZ= E AP,,* - + s+, - s-,J(29e) 

ieQ{m) 

> 0,VTO,f (29f) 

(27b), (27c)} 

where K, is the index set for uncertainty points e which are 
dynamically generated in (SP) with iterations. It should be 
noted that is the extreme point of U. Variable is 

associated with e^. The objective function in (SP) is the 
summation of non-negative slack variables ^ and s“ j, 
which evaluates the violation associated with the solution 
{x, y, z, I, P) from (MP). ^ and s“ ^ are also explained as 

un-followed uncertainties (generation or load shedding) due to 
system limitations. 


Appendix C 

Lagrangian Function for Problem (RED) 
Please check equation (30) in the next page. 


Appendix D 

Proofs for lemmas and theorems 
A. Proof of Lemma 1 

Proof Consider ^ > 0, < 0. With a small perturbation 

5 > 0 to (, we replace ^ with ^ — (5 in (RED). 
As the < 0 , then the optimal value to problem 

(RED) increases. It means that there are violations for the 
original optimal solution P^ t to problem (RED) with ^ — (5. 
Hence, the optimal solution P^ t to problem (RED) cannot be 
immunized against the uncertainty e^^ — 5. It contradicts with 
the robustness of the solution Pi^f Therefore, if i > 0, then 
> 0. Similarly, if ^ < 0, then <0. □ 


B. Proof of Lemma 2 

Proof Assume i G (z(m), according to the KKT condition 

5£(P,AP,A,a,/3,r7) _ 

d^Ptt 


at the optimal point, we have 


: - Pi + - ul)P^ = 0- (32) 

I 


and 


(SP) max min > > (sit * 

eeu is+,s-,AP)&n(e) ^ ^ ’ 




) (29a) 
(29b) 


7 ^(e):= {(s+,s-,AP) : 

E ^Pi,t = E(^™. 

i m 

-Fi < E T/.™ (pZt + ^PZ) < Fi,yi, t (29d) 


(29c) 


Then (17) holds. If nZt > 0, then ^ > 0 as 

and are non-negative. According to the 
complementary conditions for (7b) and (7d), at least one of 
(7b) and (7d) is binding. Hence, 

API = min {i.,tPr" - Pr,u Rra - ht)} 

holds. Similarly, the other equation holds when <0- 
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/:(P,AP,A,a,/3,7?) 

E E + E ^‘(E - E + E E - ktPrn+PrAktPf^ - p^A 

t i t 771 i t i 

+ EE (“*’* {^A,t - Pi,t-i - a“(i - yi,t) - PAt"'yi,tj + - Pz,t - r-f (1 - z,A - P^Am^t 


+E E (E A,mPZ - Fi) - m,t (E + ^ 0 ) + E E E - E 

t I ^ m m ^ k^K t 

+ E E E {^Up^,t+^plt - ktpm+§UktPr^ - k,t - ^pk 

kGK t i 

+ E E E {<t[^ptt - ^r(i - At)) - Qit(^pit+Rfk - 

fee/C t i ^ ^ 

E E E fe ( E k,m{kAt + AP“jf) - Fi) - rjl ( E k,„r{Pkt + AP"jf) + p,) 


fee/C t i 


(30) 


C. Proof of Theorem 1 

Proof The energy congestion cost at t is 

'y ' (km,Am,t ~ y ' ^7n,tPi,t^ (33) 

m i£Q(m) 

in 

=E E + E - iii,t - E 

m I \ keK keK / 

=E [kt + E tt] {pi - A/,7) 

I \ kGK / 

=E (+ E + VI, t + E vtt] Fi 

I \ keK keK / 

-EE(t,^fu+A^/r,‘) 

I k&K 

-Ek^fu 

I 

=E + E vtt + vi,t + E vtt] Fi 

I \ keK keK / 

-EE{tA!r:+iAiz‘) 

I k&K 

The first equality holds following the definition of net power 
injection. The second equality holds according to (8) and 
'YhmP'mt — The third equality holds following (19) and 
(20). The sign change of rji^t and ryft in third equality 
is because of the definition of power flow direction. According 
to the complementary conditions, the third term in the fourth 
equality must be zero based on the following three cases. 

1) If fii,t 7^ 0, then A/7 = 0, and rji^t = 0. 

2) If p,.t = 0 and rji^t f 0, then A/7 = 0. 


3) m,t = 0 and rji^t = 0 . 

The second term in the last equality corresponds to (21). 
The credits to FTR holders can be written as 


y ] km,t '^n,t)FTRm->.n 

(m—>-n) 

/ ^ ^ Vl,t^ \ 


(34) 


= E 

(m—>-n) 


FTR, 


I keK 

-A + E “ Vl,t) 

I 

+EEr/.n( 7*-!/7 

I fce/c 

/ (r,.„-r,,™)(i7,,, + E7t) \ 


7 


-E E 

(m—>-n) I 


k£K 


-(r,.„ - Ti,m){vi,t + E vlt) 

keK 


FTR, 


aE 


m,t + E vi,t + vi.t + E Vi,t ) Fi 
\ keK keK ) 


The first equality holds according to ( 8 ). The inequality is true 
as the amount of FTRm_).„ respects 

—Fi < E (r,,m — r,_„)FTRm_j.„ < Fi. 


according to the SFT for FTR market [18], [20], [21]. The 
right-hand-side of the inequality is the first term in the last 
equality of (33). Based on (33) and (34), the maximum 
difference between the FTR credit and the energy congestion 
cost is equal to the transmission reserve credit. That is, the 
maximum FTR underfunding is (21). 

□ 


D. Proof of Theorem 2 

Proof According to Theorem 1, the FTR underfunding value 
is ( 21 ) due to the deficiency of the energy congestion cost. 
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Therefore, we need to prove that the money collected from 
uncertainty sources can cover the FTR underfunding and 
credits to generation reserve. 

Without loss of generality, we consider the payment col¬ 
lected from uncertainty sources at time t for 

m 

m I 

-Y.Y. T/.™ {nit-ul) ( E 

i m I igp(m) 

+y:(,jf,A/,7+,*,A/,7) 

I 

=E E c‘.a7, + ^(* A/7+5‘, A/,7) 

i I 

The first equality holds according to (9). According to (7a), 
(7f), and (7g), the X)™ E/ T; ,m'riiPm t '^he second line can 
be rewritten as 

E E ( E + At) - d^,t) - E tAi 

m I I 

^EE^f-^M E ^pt,t+T.tt{T.AmPZ-Fi) 

m I ieQ(m) I m 

= EEr^™^M( E 

m I i&Q(m) I 

Em E/ t tie reformulated similarly. Hence, 
the second equality holds. The third equality holds from (9). 
Therefore, 

EE'^-.‘ = EE®m+EE®m 

m t it It 

holds. That is, the uncertainty payment covers the generation 
reserve credit and transmission reserve credit. Then, following 
the energy congestion cost shown in (33), 

m t m t 

>EE<.*^m + EE0m 

'it z t 

+ E E - <*)FTRm^n 

t m—¥n 

holds. That is, the total payments collected from loads and 
uncertainty sources can cover the total credits to energy, 
generation reserve, and FTR holders. So, the revenue adequacy 
of the proposed market clearing mechanism is guaranteed. □ 


E. Proof of Competitive Equilibrium 

Proof Pi^t and are coupled by constraints (15) 

and (16). According to (17), we can rewrite generation reserve 
credit as 


u.up^up I u.down^down 
^m,tQt,t+^m,t Qr,t 


E 


u,fe 


^Plt 


keK 


keK 


&lt)APt 


-E 

k^K 


(BUktPr^ 


p^,t)+§up^,t-ktprn\ 

t/i.i)) + Q:7t-Rf(l — ^i.t+l) / 


Substituting (35) into problem (PMP/), we can decouple Pi t 
and (Qi’’/, In fact, we also get all terms related to 

Pi t in Lagrangian C{P, X,a, f3,ri) for problem (RED). Since 
the problem (RED) is a linear programming problem, the 
saddle point Pi^t, which is the optimal solution to (RED), is 
also the optimal solution to (PMPi). Consequently, unit i is 
not inclined to deviate its output level as it can obtain the 
maximum profit by following the ISO’s dispatch instruction 
P/ t. Therefore, dispatch P^/ and price signal t^'^mt) 
constitute a competitive partial equilibrium [27]. □ 
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